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ABSTRACT 

We investigate the effect of a planet on an eccentric orbit on a two dimensional low 
mass gaseous disk. At a planet eccentricity above the planet's Hill radius divided by its 
semi-major axis, we find that the disk morphology differs from that exhibited by a disk 
containing a planet in a circular orbit. An eccentric gap is created with eccentricity 
that can exceed the planet's eccentricity and precesses with respect to the planet's 
orbit. We find that a more massive planet is required to open a gap when the planet 
is on an eccentric orbit. We attribute this behavior to spiral density waves excited at 
corotation resonances that increase the disk's eccentricity and exert a torque opposite 
in sign to that exerted by the Lindblad resonances. The reduced torque makes it more 
difficult for waves driven by the planet to overcome viscous inflow in the disk. Spectral 
energy distributions of transitional disks may reveal the presence of an eccentric planet 
if they are matched by an inner edge comprising a range of temperatures set by the 
semi-major axis and eccentricity of the hole. 

Key words: planetary systems : protoplanetary disks 



1 INTRODUCTION 

The recently discovered extra solar planets include Jo- 
vian mass planets with short period orbits and iso- 
lated planets with large eccentricities 1 . Explanations for 
the "hot Jupiters" include orb ital migration models (e.g. 
IPapaloizou and Terauernll2006t) . Explanations for the later 
inclu de eccentricity growth v ia planet /disk interactions 
(e.g. lOoldreich and Sail 1200.4 ISari and Ooldreichl l2004t 
Kiev and Dirkscn 2006), or planet /planet interaction s (e.j 



Masset and Sneilgroveil2001t iThommes and Lissauerl 1200.4 
Papaloizoul l2003l : iKlev et, alJ Eo04). Protoplanets could also 



be formed with moderate eccentricity due to interac tions 
and waves in the disk (e.g. lCresswell and Nelsonll200fift . Be- 
cause protoplanets embedded in disks might be in eccentric 
orbits we are prompted to study gap formation by an eccen- 
tric planet. 

So called "transitional disks" (systems with an accre- 
tion disc signature, where the accretion disc does not appear 
to ext end to the star) have been observed for a number of 
years <|Mar sh and Mahoney|l992: Jcnsc n~and Mathieulf 997l 
iRice et alJl2003t iBergin et all2004) . With the launch of the 
Spitzer Space Telescope, observers have been able to obtain 
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high resolution spectra of these systems (previously, only 
broadband colours were available), dramatically improving 
our understanding of these systems. In particular, Spitzer 
has established that the inner dis ks are very empty, and that 
the e dges are quite sharp (iD'Alessio et all2005HCalvet et alJ 
2005). These objects, such as CoKuTau/4, likely harbor 
massive planets residing just interior to the disk edge. In- 
teractions between the disk and planet and multiple planet 
interactions could cause this outer planet to be on an eccen- 
tric orbit. 



In this paper using 2D hydrodynamic simulations, we 
investigate the ability of a planet on an eccentric orbit to 
open a gap in a low mass viscous disk. We also investigate 
possible differences in disk morphology that are peculiar to 
planets on eccentric orbits and so would allow an observer 
to differentiate between a disk perturbed by an eccentric 
planet and one in a circular orbit. 

In section |21 our simulations are described. In section 
we discuss how we identify a gap from the numerical simula- 
tions. We discuss comparison numerical simulations of gaps 
opened by planets on circular orbits. In section[I]we discuss 
how a disk responds to a planet on an eccentric orbit. In sec- 
tion|^]we present a modified gap opening criterion, based on 
our simulations, for planets on eccentric orbits. A summary 
and discussion follows. 
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2 NUMERICAL EXPERIMENTS 

To investigate the evolution of gaps opened by a single pro- 
toplanet, we carried out a set of two dimensional hydro- 
dynamical calc ulations using the Fargo code developed by 
iMassetl i2000bl lah . Fargo is an Eulerian polar grid code with 
a staggered mesh and an artificial second order viscous pres- 
sure to stabilize the shocks. The code allows tidal interaction 
between one or more planets and a 2D non-self-gravitating 
gaseous disk It owes its name to the Fargo advection algo- 
rithm that removes the average azimu thal velocity for the 
Courant timestep limit {Massct 20 00bT) . The simulations are 
performed in the frame rotating with the planet's guiding 
centre. The outer boundary does not allow either inflow or 
outflow, so it must be located sufficiently far from the planet 
to ensure that spiral density waves are damped before they 
reach it. This is facilitated by adopting a logarithmic grid 
in radius. The grid inner boundary only allows material to 
escape so that the disk material may be accreted on to the 
primary star. While the gas disk feels a gravitational per- 
turbation from the planet, the planet itself was not allowed 
to feel the disk. This choice allows us to investigate the role 
of the planet's eccentricity in perturbing the disk without 
the additional complication of a varying planet eccentricity 
or inward migration. 

The code uses units such that G = M« = o p = 1, where 
M* is the mass of the central star, and a p is the planet's 
(initial) semi-major axis. The planet mass, M p , is described 
in terms of the ratio of the planet mass to the stellar mass, 
q/M t . The grid extends between r m i n = 0.2 and r max = 5.0, 
and a full circle in azimuth. We used 384 equally spaced cells 
in the azimuthal direction, and 200 logarithmically spaced 
cells in the radial direction. This choice allows the grid cells 
to be nearly square and so minimises truncation errors. The 
planet is initially placed at apocenter. 

We initially adopted a grid identical to that of the hy- 
drocode comparison of de Val Borro et al . , 2006 MNRAS 
in press, but found that spiral density waves were spuri- 
ously reflected by the outer boundary when the planet was 
on an eccentric orbit. We increased the maximum radius 
to that listed above and verified that our numerical experi- 
ments were not sensitive to the location of the outer bound- 
ary. 

Our disk viscosity is constant over the entire disc. We 
typically choose v — 10~ in the system of units outlined 
above (although we have performed runs with other values). 
Once the viscosity is set, we may calculate the Reynolds 
number of the disc: 

In our units, TZ = v^ 1 at the radius of the planet's orbit 
(note that TZ is a function of radius). 

We use a constant aspect ratio disc, with h/r — 0.05, 
where h is the density scale height. The local sound speed 
is set from this. For our initial conditions, we took a flat 
surface density profile (E(r) = Eo). Since the planet was 
not allowed to migrate, Eo was only important numerically; 
we verified that it was sufficiently small to avoid spurious 
numerical effects. 
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Figure 1. The azimuthally averaged surface density of the disk 
after 500 orbits for planet mass ratios 10 — 4 < q < 10~ 3 . All runs 
had the planet on a fixed circular orbits, and the disc viscosity 
was v = 10 -5 

3 GAP OPENING IN THE CIRCULAR CASE 

Planets above a certain mass are expected to open a gap in 
their natal disks. In this section, we shall discuss the for- 
mation of gaps by planets on circular orbits. Although this 
has been done by numerous authors, it is important that we 
first characterise the behaviour of Fargo for planets in a 
circular orbit, before we proceed to the eccentric case. 

Previous studies have shown theoretically and numeri- 
cally that a gap should be opened in a disk when the planet 
mass ratio exceeds 

q > 40ft -1 (2) 

<Lin and Papaloizoul Il993 iBrvden et al]ll999h . This con- 
dition can be derived by balancing the viscous torques, 
T„ = 3TvT,ur 2 Q,, against those from the spiral density waves 
the planet excites at its Lindblad resonances. 

In Figure Q we show the azimuthally averaged surface 
density after 500 orbits for planet mass ratios ranging from 
q = 10 -4 to 10 -3 on fixed circular orbits. These runs all 
have v = 10 -5 for their viscosity. We see that the depth 
of the gap varies smoothly with q, so determining whether 
a gap has formed is somewhat arbitrary. Evaluating equa- 
tion |21 we find that q — 4 x 10~ 4 is the expected threshold 
mass ratio for gap formation for this set of runs. Figure 
shows that mass ratios larger than this have density con- 
trasts E gap /Eo < 0.2 at the bottom of their induced gaps. 

Having tested the behaviour of the code as q varies, we 
now turn our attention to the variation with viscosity. Fig- 
ure [5] is similar to Figure except the planet mass ratio is 
fixed at q — 10 -3 , and we varied the disk viscosity. Equa- 
tion [5] predicts that for a protoplanet with q = 10 -3 , a disk 
with Reynolds number at the planet above TZ = 8 x 10 4 
should cause a gap to open in the disk. This corresponds 
to a viscosity of v < 2.5 x 10 -5 . The runs satisfying this 
criterion all show E gap /Eo < 0.2 for their gaps. 

Based on the surface density profiles shown in Figures Q 
and|5J we adopt the following criterion to identify formation 
of a gap in our numerical experiments: the surface density 
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Figure 2. The azimuthally averaged surface density of the disk 
as the viscosity varied between v = 10 — 5 and 7x 10 — 5 . The deeper 
gaps correspond to the lower viscosity disks. All simulations are 
shown after 500 orbits, for a q = 10~ 3 planet on a circular orbit 




Figure 3. The azimuthally averaged surface density of the disk 
after 500 orbits for planet eccentricities between < e p < 0.12. 
The planets had q = 6x 10 — 4 and the disc viscosity was v = 10~ 5 



in the gap must be less than 20% of the unperturbed surface 
density. We have chosen this criterion so that it is consis- 
tent with equation [5] We adopt this criterion during our 
discussion of the effects of eccentric planets in subsequent 
sections. However, as Figures0and|5]show, the precise def- 
inition of 'gap' is arbitrary, since the variation with q and v 
is relatively smooth. 



4 DISK RESPONSE TO AN ECCENTRIC 
PLANET 

In this section we discuss how a disk responds to a planet 
on an eccentric orbit. We concentrate on planets with mass 
ratio q = 6 x 10 -4 and 10 -3 in disks with viscosity v — 10 -5 . 



4.1 Azimuthally averaged density profiles 

In Figure [3] we show how the surface density profiles vary 
with planetary eccentricity, e p for ag = 6xl0~ 4 planet. 
We see that the shape of the gap is similar to the circular 
case for planet eccentricity below e < 0.06. The planet's Hill 
radius is given by 




Figure 4. Similar to Figure [3] but the planet is mass ratio 
q = 10 -3 and the eccentricity range is < e p < 0.2 



(3) 



which measures the distance over which the planet's gravity 
dominates (the Hill sphere is the size of the Roche lobe in 
the low q limit). In our units, we note that rnm = 0.058 for 
these runs. Above this eccentricity, the gap as seen in the az- 
imuthally averaged density profiles, is shallower and wider. 
By the 20% criterion discussed in the previous section, the 
planet no longer induces a gap for e p > 0.1. 

Figure0]is similar to Figure[3]but for a q — 10~ 3 planet. 
At eccentricities of e p < 0.065, the shape of the gap is nearly 
identical to that of the circular case. Above this threshold, 
the gap becomes shallower and wider, and increasingly de- 
viates from the circular case. We note that the planet's Hill 



radius is 0.069a p . This suggests that the planet's eccentric- 
ity has little effect on the density profile unless the planet's 
orbit takes it outside the sum of its semi-major axis and Hill 
radius. As spiral density waves are driven from resonances 
located at radii all the way up to the planet's Hill radius, 
this limiting eccentricity is not unexpected. 

Figures |^1 and 0] show that as planet eccentricity in- 
creases, the azimuthally averaged gap density profile be- 
comes wider and shallower. It becomes harder for the pro- 
toplanet to open up a deep gap in the disk. Because the 
density in the gap is lower at higher planet eccentricity, a 
more massive planet is required to open a gap in the disk. We 
shall discuss possible explanations for this behaviour after 
we describe the disk morphology seen in these simulations. 
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4.2 Eccentric gaps 

In Figure |K| we show the 2D surface density of the disc 
after 250 and 500 orbits for planet eccentricity e p = 0.2 
and mass ratio q = 10~ 3 . Examining the surface density 
as a function of time, we see the planet crossing between 
the inner and outer gap edges, producing one armed spirals 
as it encounters each edge. By t — 500 orbits, the planet 
has formed an eccentric gap. The eccentricity of the gap 
increases to a maximum of about twice that of the planet at 
this time. The gap also slowly precesses. At t = 500 orbits 
the gap's apocentre is about 180° from that of the planet. 

Because the gap precesses with respect to the planet's 
orbit, the planet approaches the disk edge at different lon- 
gitudes. Close approaches allow the planet to clear out ad- 
ditional material from the disk edge. On longer timescales, 
t > 1000 orbits, the gap slowly widens and eventually circu- 
larizes. However during the long circularization period ac- 
cretion onto the planet could exceed that of a similar mass 
planet on a circular orbit. In the case of a planet on a cir- 
cular orbit, the formation of a gap can significantly reduc e 
the accre tion rate onto the plane t flD'Angelo et al.l l2003). 
Recently, iKlev and Dirksenl (|200(t) looked at the excitation 
of disc eccentricity by massive planets (up to q = 5 x 10 -3 ), 
and suggested that the induced eccentricity could aid accre- 
tion by the planet. Although we did not permit accretion in 
our numerical experiments, we would expect behaviour of 
this nature. 

Figure 0] showed the azimuthally averaged density pro- 
file was shallower for gaps opened by eccentric planets. The 
smoothness or shallow gradient of the azimuthally averaged 
profile in the gap is due to the eccentricity of the gap edge, 
rather than a smooth change in the disk density with radius. 
The gap itself has sharp edges. The apparent shallow profile 
in the azimuthally averaged profiles is caused by averaging 
the density profile at different longitudes. 

Figure |S| plots the eccentricity evolution of the gap for 
three values of e p . To determine the eccentricity of the gap, 
we selected cells based on the E/Eo < 0.2 criterion intro- 
duced in section |21 We then computed the cell eccentricity 
by treating each as a free particle, and took the mean of 
these. Because of the rather arbitrary gap criterion used, 
the time evolution shows some large oscillations. 2 The cir- 
cular (e p = 0) case shows a brief initial spike of eccentricity 
in the gap, before decaying to a low value. For the two runs 
with an eccentric planet, we see a steady growth in gap ec- 
centricity, approaching 2e p after 1000 orbits. 



5 GAP OPENING AS A FUNCTION OF 
PLANET MASS AND ECCENTRICITY 

We performed a series of runs, all with v — 10~ 5 , and with 
different planet masses and eccentricities. For each run we 
measured the depth of the density deficit near the planet and 
labeled the simulation as "gap opening" if the density near 
the planet was lower than 20% of the unperturbed value. 
Figure |7| shows the result of these runs. We distinguish be- 
tween three behaviours: no gap opened (by the 20% criterion 

2 We could produce a weighted mean, but the appropriate weight- 
ing is not obvious, so we prefer to show the simple average 
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Figure 5. Disc surface densities after a) 250 and b) 500 orbits 
for a q = 10~ 3 planet on a orbit with eccentricy fixed at e p = 0.2. 
The gas viscosity was v = 10~ 5 

introduced above) ; a gap formed which is significantly differ- 
ent to the circular case; and a gap formed which is essentially 
identical to the circular case. 

From Figure |7| we see that simulations similar to the 
circular case are found at the bottom of the plot. The di- 
vision between density profiles resembling the circular case 
and those differing is not a strong function of eccentricity. 
The border between the two regions is well described by 

^threshold ~ (^) 

This is essentially equation [3] and compares the variation in 
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Figure 6. The eccentricity of the gap opened in the disk as a 
function of time for a q = 10 — 3 planet, for e p = 0, 0.1 & 0.2, 
and a disk viscosity of v = 10~ 5 . We computed the eccentricity 
by locating all the cells satisfying E/Eo < 0.2, treating each as a 
free particle, and taking the mean of the resulting eccentricities 
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Figure 7. Gap formation as a function of planet mass ratio and 
eccentricity. For these runs the disk viscosity was v = 10~ 5 so a 
gap is expected for q > 4 X 10 — 4 . The shaded square points refer 
to runs in which a gap was opened in the disk. The empty square 
points represent runs that closely followed their respective circu- 
lar cases when a gap was formed. The X indicates a run lacking 
a gap. The dotted line shows the planet's Hill radius divided by 
the planet's semi-major axis. Simulations below this line closely 
resembled their circular orbit counterparts. At large planet eccen- 
tricities, the planet could not open a gap. We have empirically fit 
a function, shown as a dashed line to mark the division between 
gap opening simulations and those lacking a gap. This function is 
in the form (equation|SJ expected if co-rotation resonances reduce 
the total torque on the disk due to density waves driven by the 
planet 



is not opened in the disk even when equation [5] would pre- 
dict one. We now discuss possible reasons for this. We ex- 
pect a gap is opened when torque due to damping of den- 
sity waves exceeds that due to viscous flow. We first con- 
sider the accretion rate throug h a viscous and eccentric disk. 
ISver and Clarke! l|l992t E$93) showed that the torque due 
to a viscous accretion varies with eccentricity according to 
T v oc (1 — e 2 ) (see their equation 2.16). However this reduces 
the viscous torque rather than increasing it, and would make 
it easier for an eccentric planet to open a gap rather than 
harder as we see from Figure |7| 

The solution lies in the other reseonances present in 
the system. Because of perturbations from the planet, spi- 
ral density waves are excited at both Lindblad and co- 
rotation resonances. Lindblad resonances are expected to 
increase the eccentricity of the planet, a nd hence would de- 
creas e the eccentricity of the disk (e.g. ISari and Goldreicbl 
12004) . Co-rotation resonances, however, when in the pres- 
ence of a density gradient, damp the eccentricity of the 
planet and hence are expected to increase the eccentric- 
ity of the disk. First order terms in the potential, propor- 
tional to the planet eccentricity, can drive co-rotation reso- 
nances into the disk at the same location as the traditional 
Lindblad resonances. Notably the corotation resonances are 
also expected to have torque with the opposite sign as the 
Lindblad resonances. Waves driven at these resonances and 
damped in the disk do not push the disk away from the 
planet, but rather pull the disk toward the planet. Conse- 
quently the corotation resonances reduce the torque on the 
disk due to spiral density waves driven by the planet. The 
torque from these waves Ten o c ^7? 2e p (from equation 14 of 
iGoldreich and Tremainelll980l) We therefore expect a mod- 
ified gap opening criterion of the following form 
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where \ depends on the density gradient and the ratio be- 
tween the sum of torques from the corotation resonances and 
the sum of torques from the Lindblad resonances. Unfortu- 
nately \ is dependent on the density profile and is not easy 
to calculate. The parameter \ may be large since a sharp 
gradient in the density could cause a large torque even if the 
gap is not deep. Since the shape of the gap is certain to be a 
function of viscosity (cf Figure we expect x to vary with 
viscosity too. However, we did not vary viscosity in this set 
of runs. 

We have fit a function in the form equation |5] consistent 
with the gap/no-gap line seen in Figure |7| The gap/no-gap 
curve is consistent with x ~ 912. This suggests that the 
corotation resonances are primarily important in regions of 
strong density gradients. Only a planet massive enough to 
cause a moderate density gradient in the disk would be able 
to increase the disk eccentricity. 



peri- and apastron to the planet's Hill sphere. This suggests 
that the planet's eccentricity has little effect on the den- 
sity profile unless the planet's orbit takes it outside the sum 
of its semi-major axis and Hill radius. Spiral density waves 
are driven by resonances up to the edge of the planet's Hill 
sphere (this assumption is used to derive equation^, so it is 
not unsurprising that once the planet strays over these reso- 
nances (effectively destroying them) the behaviour changes. 
Figure |7| shows that at high planet eccentricity a gap 



6 DISCUSSION 

We have seen that an eccentric planet can open an eccen- 
tric gap in a disc. The eccentricity should be an observable 
quantity. Although we see eccentricit y in images of dis cs - 
the eccentric ring around Fomalhaut llKalas et alj|2005l) for 
example - we would expect eccentric holes of the nature de- 
scribed here to only be a few AU from their host stars. Such 
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holes would be too small to be imaged by current observato- 
ries. However, a strongly eccentric gap should leave a mark 
on the spectral energy distribution (SED) of a disk. The sur- 
face temperature of a disk illuminated by the star is set by 
the radial distance T oc r -1 ^ 2 . Consequently the eccentric 
edge of a disk would correspond to a range of dust tempera- 
tures in the edge. The spectral energy distribution of CoKu- 
Tau/4's disk edge w as consistent with a si ngle temperature 
optically thick wall (D'Alcssi o et, al]l2005r) . However an ec- 
centric hole would produce a spectral energy distribution 
consistent with a range of dust edge temperatures. Such a 
spectrum could either be interpreted as due to a smooth ra- 
dial density gradient or an eccentric hole. If the disk were 
optically thick then the second possibility would be more 
likely than the first. We have found that the disk edge can 
be more highly eccentric than the planet, consequently the 
variation in edge temperature around an eccentric disk edge 
may be detectable even though the dust temperature is only 
proportional to the inverse of the square root of the radius. 

How might such a system form? In our r uns, the planet 
was on a fixed orbit. IPapaloizou et alJ (1200 ll) note that the 
planet and disc should exchange eccentricity, with the ex- 
change being most efficient when the planet and disc are of 
comparable mass. A massive disk would also be exp ected 
to damp the planet's eccentricity l)Artvmowiczl I1993T) . We 
are therefore lead to two possible situations in which an 
eccentric gap might arise. The first is the case of an ec- 
centric planet in an old, dissipating disk. In this case, we 
expect the planet to be more massive than the disk, and 
hence little affected by it. How the planet achieved its high 
eccentricity is then an open question. The second scenario 
requires a second, inner planet. If locked in resonance with 
the outer planet, it could 'resupply' eccentricity to the outer 
planet, even in the face of damping by the outer disk. We 
might even expect the outer planet to have been pushed 
into the resonance by conventional planet-disk migration. 
This is so mewhat analagou s to the case of GJ 876, as dis- 
cussed bv lKlev et alJ i2005l) . However, we would expect the 
inner planet to be more massive, and hence less affected by 
the resonant interaction. This is not the case in GJ 876, 
where the inner planet is believed to b e lower ma s s, and 
on a substantially more eccentric orbit. iKlev et al.1 (120051) 
required rapid dissipation of their gas disk, to prevent the 
inner planet becoming too eccentric. 



7 SUMMARY 

In this paper, we have examined the effect of an eccentric 
planet on a gas disk. From our numerical experiments, we 
propose that the conventional gap-opening criterion should 
be modified to become 

g(l-912ep) = 40ft" 1 (6) 

The form of equation |S| is justified by consideration of coro- 
tation torques, which are expected to act to oppose the 
Lindblad torques. However, the extra coefficient had to be 
determined numerically, due to the sensitivity of the coro- 
tation torques to the exact gap profile. We expect its value 
to depend on the assumed disc viscosity, but we have not 
parameterised this variation. 



We have found that the gap eccentricity can substan- 
tially exceed the planet's eccentricity. In the early stages of 
evolution, as the planet crosses between the inner and outer 
discs, it induces one armed spirals in each. Coupled with the 
precession of the gap, we expect this to lead to accretion by 
the planet, despite the presence of a gap. 

If planet's eccentricity is sufficiently low, then the gap 
formed is almost identical to the circular case. We have 
found the threshold eccentricity to be well fit by 



which compares the peri- and apocentre of the planet to 
its Hill sphere. The form of equation |7| is compatible with 
the standard theory of gap formation. This theory assumes 
that resonances up to the edge of the Hill sphere contribute 
to opening the gap (within the Hill sphere, the resonances 
do not exist, since there material orbits the planet, not the 
star). Once e p is high enough to let the planet move outside 
its notional Hill sphere, these resonances will be destroyed, 
and we would expect the morphology of the gap to change. 

We would expect an eccentric gap to be detectable in 
a disk SED. However, without corroborating data, it might 
be difficult to disentangle eccentricity from an azimuthally 
symmetric radial density gradient. 

Future work is needed, allowing the planet to feel the 
disk. This will allow us to determine the maximum disk 
mass which does not rapidly circularize the planet. We have 
also suggested that a high planetary eccentricity might be 
maintained by a second, massive, inner planet, even in the 
presence of strong damping by a higher mass disk. Further 
calculations are needed to examine this possibility. 
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